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The selection of inflow values at boundaries far upstream of an aircraft is considered, for one- and two- 
equation turbulence models. Inflow values are distinguished from the ambient values near the aircraft, which 
may be much smaller. Ambient values should be selected first, and inflow values that will lead to them after 
the decay second; this is not always possible, especially for the time scale. The two-equation decay during the 
approach to the aircraft is shown; often, the time scale has been set too short for this decay to be calculated ac- 
curately on typical grids. A simple remedy for both issues is to impose floor values for the turbulence variables, 
outside the viscous sublayer, and it is argued that overriding the equations in this manner is physically justified. 
Selecting laminar ambient values is easy, if the boundary layers are to be tripped, but a more common practice 
is to seek ambient values that will cause immediate transition in boundary layers. This opens up a wide range 
of values, and selection criteria are discussed. The turbulent Reynolds number, or ratio of eddy viscosity to 
laminar viscosity has a huge dynamic range that makes it unwieldy; it has been widely mis-used, particularly 
by codes that set upper limits on it. The value of turbulent kinetic energy in a wind tunnel or the atmosphere 
is also of dubious value as an input to the model. Concretely, the ambient eddy viscosity must be small enough 
to preserve potential cores in small geometry features, such as flap gaps. The ambient frequency scale should 
also be small enough, compared with shear rates in the boundary layer. Specific values are recommended and 
demonstrated for airfoil flows. 


I. Introduction 

Turbulence models based on transport equations are now omnipresent, but the CFD community’s knowledge of 
appropriate boundary conditions for the turbulence variables at distant boundaries is incomplete. There is much room 
for clarity and helpfulness in some of the original turbulence-model papers, and in the manuals and Graphical User 
Interfaces of widely-used codes. Solid physical interpretations of the values taken by the primary variables and the 
eddy viscosity derived from them are needed. A contributing factor is the large dynamic range of the variables; for 
instance, the eddy viscosity may logically be compared with the molecular viscosity or with the product of velocity and 
length of the aircraft, and these two numbers may differ by a factor of 10 8 , the Reynolds number. The choice between 
these two scales is a source of confusion, and interferes with the establishment of guidelines. This situation creates 
pitfalls, among them failures to produce turbulence where intended by the user, 1 and excessive eddy-viscosity values 
in regions which have not produced turbulence and should remain irrotational. These pitfalls, which are preventable, 
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may well go undetected when the models are treated as black boxes. Our purpose here is to recommend effective 
prevention measures. “Effective” refers to boundary-condition practices which do not add much to the users’ burden, 
and ensure a low sensitivity of the results to trivial parameters. For instance, it would be unacceptable to get different 
results from a computation with a domain 50 chords in radius, and another domain 60 chords in radius. Progress in 
science and engineering depends on repeatability, and on building an experience base that is genuine. 

It may help to start from the “ideal situation” a trustful CFD user would expect. The incoming stream has well- 
defined turbulence kinetic energy and dissipation, k and e, and eddy viscosity derived from them, v f = C' ; ,fc 2 /c (for 
a physical reason and not only dimensional analysis). These quantities correctly respond to velocity gradients while 
approaching the aircraft, and couple meaningfully with the boundary-layer turbulence. It is simply a matter of knowing 
k and e far upstream, preferably from measurements. These values are “realistic.” 

This vision fails, primarily because describing all of turbulence with two numbers is insufficient. Here, we disre- 
gard internal flows and in particular gas-turbine applications, for which the coupling with ambient turbulence is much 
stronger, and address only external aircraft flows. In the atmospheric boundary layer, modeled with Reynolds-averaged 
Navier-Stokes (RANS) equations, the eddy viscosity on a windy day can well take values such as 50m 2 /s, and the 
chord Reynolds number of a large airplane based on this eddy viscosity would be near 5. This is untenable, of course. 
The reason is that the eddies which may couple with the boundary layer have sizes of the order of Icto, rather than 
100m, so that only a tiny fraction of the “atmospheric” k is relevant (the analogous estimates in a wind tunnel are 
not as immensely different, but the physical consequences are the same). This relevant fraction is deduced from the 
atmospheric dissipation rate, much better than from the atmospheric k, and it is not sharply defined. Moreover, it has 
no reason to obey the decay equations used to calibrate two-equation models in isotropic turbulence. In reality, the 
kinetic energy relevant to the aircraft flow varies very little over the size of the CFD domain. This fact is in strong 
conflict with common practices in CFD, since inflow values of e/k near U^/c are routinely used, where is the 
speed of the aircraft and c the wing chord. Over an approach to the aircraft of length 50c, k then drops by a factor of 
about 75, if the solution is accurate according to the decay equations. This is obviously unrealistic behavior, accepted 
in an attempt to obtain the desired behavior in the boundary layers. Furthermore, only the eddy viscosity has a strong 
influence on the flow field, at least at moderate Mach numbers; k and £ separately do not. In fact, most CFD users 
are only interested in whether the solution is laminar or turbulent, and of those, the vast majority have no interest 
in laminar regions. Fully-Turbulent (FT) behavior is safer and gives better convergence, and can be viewed as the 
standard mode of use of RANS models today. 

In addition to being unrealistic, the precipitous decay associated with such inflow values is not calculated accu- 
rately. A typical grid spacing near the inflow boundary is larger than c, which is too large to resolve a decay at a 
rate 1/c. This is observed in otherwise fully normal CFD solutions; the decay is grossly under-estimated and grid- 
dependent, as will be shown below. The boundary-layer turbulence responds only to what we will call the ambient 
values of the turbulence variables: those in the vicinity of the aircraft. Thus, for example, grid-dependent ambient 
values can have an influence on the location in the boundary layer where the model trips to turbulence. Because small 
variations in transition location do not typically have a glaring impact on the global results, this problem is largely out 
of sight for the user. However, the situation is clearly problematic from the point of view of consistency. 

Most CFD codes handle the issue of freestream turbulence differently. There does not appear to be a consistent 
standard practice. Many codes assign default values and allow the user to override if desired. However, there is usually 
very little guidance given as to how different choices affect the decay rates and the resulting levels that are actually 
seen near the body. Some codes (such as the CFD++ code 2 ) allow the user to disable the decay in some regions, 
which are manually defined with controls such as a; and ;;:, nax . This manner of control is not practical for real- 
life geometries. Other groups prevent the decay by imposing “floor values:” the turbulence variables are limited by 
non-zero minimum values. This limiting may conflict with the wall boundary condition for k, but this can be resolved 
using the SST model’s F\ function, 3 for instance. More important is the question of whether violating the celebrated 
k-e equations by imposing floor values can be justified at all; this is taken up below. 

In the remainder of the paper, most of the focus will be on the two-equation k-u> SST model 3 and the one-equation 
model of Spalart & Allmaras (S-A), 4 but theoretical aspects of the k-e model will also be considered. The equations 
for free decay are presented first, leading to “reverse-engineering” inflow values from ambient values; constrained 
decay is then examined. Then, the valid criteria that exist to pick ambient and/or inflow values are discussed. This is 
followed by examples which support the recommendations made, and a summary. 
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II. Decay Equations and Control 


A. Free Decay 

Free decay rates of turbulence quantities can be very large and “non-physical” when commonly-used inflow values are 
specified. This can be seen from the solution to the k-e equations in the approaching flow field. A typical form of the 
complete k-e equations can be written: 
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with v t = C^k 2 /e, and = 0.09. In the approach to the aircraft, only the destruction terms are active, so that 
Uk x = —e and Ue x = -C e2 e 2 /k to a very good approximation. The solution is simple and linear for k/e: 
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where the subscript FS means “free-stream”, or “inflow”, and x is the distance from the inflow to the field point. 
Control over the time scale k/e through the inflow value is difficult; nothing can prevent it from growing. Typically, 
C s2 = 1.92. 

The solutions for k, £, and v t are then: 
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It is seen that v t decays much more slowly than k and e: the power is about —0.09 versus —1.1 and —2.1, respectively. 
Numerical tests have confirmed this decay, but have also shown that it is calculated accurately from the onset only 
if ( e/k ) is small enough compared with (Too /Ax, where Ax is the streamwise grid spacing. When Ax is too large 
to support accurate computations in the far field, the decay simply tends to be delayed, and during the approach, 
both (e/k) and Ur^/ Ax evolve in the direction needed to improve their ratio. This situation remains aesthetically 
undesirable, and likely to cause jumps during grid sequencing, among other things. It could also disturb the order of 
accuracy of the solver. 

The freestream decay of u> (for the k- to model) can be found in a similar fashion. For example, the k-to SST model 
is given by; 
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where v t = aifc/max(aiw, klF 2 ), a\ = 0.31, f 1 is the magnitude of vorticity, and Fi and F 2 are blending functions. 
Here we use the standard convention that u> = e/(C^k). The decay of o o in the freestream is: 


uj = u FS (l + w fs /3x/U). 1 (9) 

Because /3 is equivalent to C^(C £ 2 — 1) (with C p = /3* = 0.09), Eq. (9) is equivalent to Eq. (3), which can also be 
written as 1/w = 1 /ojfs + C fl (C e 2 — 1 )x/U. 

Equations (4 - 6) and (9) are available to produce the inflow values needed, once ambient values have been chosen, 
or else to indicate that a set of ambient values is unreachable, with (3) expressing the dominant constraint. 

Clearly, it is seen from these equations that the chosen freestream values of the turbulence quantities have a strong 
influence on the rate of decay of the turbulence values. Table 1 shows examples of typical freestream levels (for 
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Table 1. Examples of effect of freestream levels on turbulence decay after x/ L = 50 (for Re f = 10 million) 
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Figure 1. Evolution of ujLjU ahead of NACA 0012 airfoil, with Re c = 10 7 . 
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ReL = 10 x 10 6 ), along with their decayed levels according to Eqs. (4 - 6) and (9) at a distance of x/L = 50 
downstream. The freestream turbulence intensity, in percent, is related to turbulent kinetic energy by the relation 
Tu = 100\/2fc/(3[/ 2 ). The highest value u>L/U can reach is 0.2415; many of the inflow-value sets chosen lead to 
0.24, in other words, u>fs was set so high it lost control of the decayed value. The dramatic decay of k and e is evident 
in some of these examples. Taking Tu and ft,FS as the parameters easiest to relate to and varying them independently, 
this table demonstrates that the higher the freestream turbulence intensity (Tu), the more rapid the decay; and the lower 
the freestream eddy viscosity level (z/ t ), the more rapid the decay. In other words, to achieve “reasonable” decay rates, 
one is forced to use either very low levels of freestream Tu or very high levels of freestream v t . 

As mentioned earlier, for large decay rates, the coarseness of the grid far from the body may make the solver 
unable to support an accurate computation and the decay can be grossly underestimated. The behavior of ui for 
two different freestream turbulence levels is shown in Fig. 1. In particular, it is shown that for the higher freestream 
decay conditions, the computed results display a very noticeable grid dependence, only slowly approaching the correct 
answer as the grid is refined (upper curves). Using significantly lower values for freestream u>, however, the decay rate 
is dramatically reduced, and CFD predicts the levels with greater precision on all three grids (lower curves). 

The equivalent situation for the S-A model is most simple; in the absence of significant destruction term (which 
is discussed below) and of vorticity (which sets production), the transported variable v and the eddy viscosity v f 
are uniform. The objection sometimes made to the absence of a destruction term, which may appear unphysical, is 
here a definite advantage. A trivial issue which nevertheless causes confusion is that some user interfaces do not 
make it clear whether v or v t is being set; unfortunately, this is pivotal when the values are of the same order as the 
molecular viscosity v. An obliquely related issue is that several groups have used for inflow the value /> = 1.34, 
which gives v t /v = 0.009; this comes from a historical detail of the Baldwin-Barth model, 5 using the round number 
Ret = k 2 jve = 0.1. Unfortunately, this value is not very far from the boundary vjv = yjlog^ctz) / CtA between those 
solutions that are attracted by 0 once the production term becomes active, due to the f t 2 function in the model, and 
those that are not. 1 This boundary is at vjv = 0.60 in the more common Version la, 4 and at vf v = 0.22 in Version 1 
(AIAA paper 92-0439). The destruction term can steer v into the basin of attraction, with unpredictable results. Values 
of v/v in the 3-to-5 range are much more reliable when fully-turbulent behavior is desired. 

B. Constrained Decay 

The concept of constrained decay - or “floor values” - is simple: k and e (or ui) are not allowed to drop below some 
pre-set values. This method is used in the NASA PAB3D code, 6 and was suggested to us by Dr. Abdol-Hamid (2006). 
For the S-A model, this technique is not necessary. 

As mentioned earlier, the conflict with the wall boundary condition k = 0 can be removed by conditioning the 
limiters on F\ from the SST model, or a similar quantity in other models that detects the interior of the boundary layer. 
This conditioning should also prevent potential problems for cases with boundary layer relaminarization. A secondary 
issue is that preventing the free decay with an over-ride can create the appearance of non-convergence of the solution, 
since the residuals will remain non-zero; again, it is remediable and only a reporting issue. 

The objection against disabling the classical k-e equations is deeper, but the discussion above regarding the ac- 
tual behavior of that part of the turbulence that is of small scale compared with the atmospheric scales or even the 
wind-tunnel size is founded. This small-scale component of a larger turbulence field simply does not decay the way 
turbulence that has all its energy in the small scales would. Essentially, it has a supply of energy in the Kolmogorov 
cascade, and does not decay. A separate argument is that the physical merit of any influence the ambient values have 
inside the boundary layers via the model’s diffusion terms is very debatable. In fact, this influence has at times been 
grievously troublesome, in the case of the pure k-u model, and of the Baldwin-Barth model. See Menter. 7,8 

A milder version of this “re-interpretation” of k and e as not representing the totality of the turbulence occurs in the 
logarithmic layer. The two equations demand a plateau with k + = 1/ \JC t , , yet the true TKE found from experiments 
and DNS is far from uniform, leading to the view that the modeled k represents only the “active” energy, which has 
length scales of the order of the wall distance and contributes to the shear stress. In that sense, it is also a high-pass 
filtered version of the complete k. The description of the k-co model by Wilcox 9 also hints at such a view; see also v' 2 
in the v 2 — / model. 10 

At first sight, setting floor values adds to the user’s burden; but in fact, these levels can be identical to the inflow 
values, and they become the ambient values. Thus, the situation is markedly simplified, and this option is practical. 
We can now turn to the choice of these ambient values. 
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III. Ambient Values 


Turbulence modeling is a pragmatic endeavor, and we are seeking “effective” values in the sense outlined above, 
without reference to “true” measurable levels of k and s, let alone v t . We also consider the common situation in which 
the ambient values are picked to trigger turbulence in all boundary layers without local tripping, except perhaps at 
very low Reynolds numbers. Before making our recommendations, we review the possible bases for picking ambient 
values. 

A. Turbulent Reynolds Number 

It is natural to compare the eddy viscosity v t to the molecular viscosity /./, but serious mis-conceptions have ensued. 
First, the idea that outside boundary layers it must satisfy v t /v 1, “since there is no turbulence,” is erroneous as will 
be shown. Values much higher than this remain “effective,” when the flow Reynolds number is high enough, unless 
of course laminar boundary layers are desired. Second, objections to very high values of v t /v are also erroneous, 
and have caused trouble in at least one code. The ratio v t /v had been limited to 10 5 , at first sight a very large 
number, in response to some convergence troubles; however, 10 5 * is routinely exceeded in jets at full-scale Aerospace 
Reynolds numbers. Any such hard-wired limiters should be removed from all codes. They create a Reynolds-number 
dependence that runs completely against turbulence theory in free shear flows. This also applies to limiters offered to 
users in some codes. In some situations, RANS models err by having too high an eddy viscosity; examples include 
free vortices, for most models, and fully-developed round jets, for the S-A model. A more accurate solution can then 
be obtained by manually lowering the eddy viscosity, and some users experiment, say with a limit such as v /v < 10 4 . 
The danger is that in free shear flows, this limiter should scale proportionally to the flow Reynolds number, and the 
user’s manuals give no guidance on this. Massive inaccuracies can again result at full-size Reynolds numbers. 

The area in which u t /v has control and is meaningfully compared with unity concerns transition in the S-A model. 
Ambient values of v/v well below 1, such as 0.1, will cause the model to remain dormant. Values well above 1, such 
as 3, will cause it to grow turbulent in all boundary layers at reasonably high Reynolds numbers. This is by design. 
There is no reason to pick values much above 3, even when fully-turbulent behavior is desired. This will ensure that 
ambient values do not influence the interior of the boundary layer, except for steering it towards laminar or turbulent 
levels, which is binary. 

B. Boundary-Layer Levels, and Device Reynolds Numbers 

If the ambient values, which have been seen to be quite arbitrary, are not to influence the boundary layer through 
diffusion, these values need to be decidedly lower than values in the boundary-layer’s outer region. These are of 
the order of 10 _3 (7 2 for k, where U is a typical edge velocity, and 10 ~ 4 U 3 /S for e, where 5 is the boundary-layer 
thickness; finally, u> is simply of the order of U/8, and v t of the order of 10 ’'Lid. It must be remembered that v t itself 
is not diffusing in two-equation models, and its ambient value is not directly relevant: the k-e equations tend to draw 
v t down near the boundary-layer edge. 1 1 However, if floor values constrain the primary variables, they also constrain 
V t . 

An example of the corruption of the solution in the boundary layer by high ambient turbulence is shown in Fig. 2. 
Here, the floor value for k/U 2 is 1.5 x 10 -6 , and the floor value for uiL/U is varied. When it is very high (over 1000), 
the natural decay of lo is interrupted only 1/3 of the way up the boundary layer, and the velocity profile gets damaged 
as a result. This illustrates a danger of the practice of imposing floor values: recall that with free decay, luL/U could 
not exceed 0.24. This figure, incidentally, confirms the estimate of U/8 for the outer-layer levels of a o. Also note the 
large dynamic range of u>, and therefore e, within a boundary layer; wall values would not be relevant in this context. 

Another very concrete consideration in aerospace applications is the gap regions between different elements of a 
wing. These are normally designed to allow a “potential core” between the two boundary layers. It is all too easy 
to set the ambient v t so high that this potential core is smeared, and to be ignorant of this fact, deeply hidden in the 
flow field. Yet, this will make CFD inaccurate, and sensitive to ambient values. The gap d is often of the order of 
1% of the chord, and the gap/eddy-viscosity Reynolds number dU /v t must be in excess of about 100, and therefore 
the chord/eddy-viscosity Reynolds number must be in excess of 10 4 . Even single-element wings have small features 
for which the potential flow should not be corrupted; the leading-edge radius is a few percent of the chord. Again, 
only extensive examinations of the flow field would exhibit such violations, so it is highly desirable to prevent them 
from ever occurring in the first place by good standard practices. Normally, the constraint for the eddy viscosity to be 
dominated by the boundary-layer eddy viscosity v t < 10 :i Ud is stronger than the gap constraint v t <C Ud, but both 
should be kept in mind. 
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Figure 2. Effect of ambient ui on the boundary layer of a NACA 0012 airfoil, with Re c = 10 1 


Figure 3 shows a solution of a multi-element airfoil at Re c = 9 x 10 6 , with the S-A model and different ambient 
turbulence levels. As v/v is increased, the contamination in the solution gets progressively worse; note, however, that 
the solution degradation only starts to become noticeable at very high levels, greater than 500 or so (the low ambient 
value of F/i/ = 5 is close to the recommended value of 3, and could be considered as “standard”). In other words, the 
acceptable range for v/v is very wide. However, note that this range shrinks with lower Reynolds numbers. 

Figure 4 displays the eddy viscosity, comparing a standard ambient value of vps/ v = 5 with a high value of 
vfs/v = 5000. Using the ambient turbulent value of 5 makes v/v visible on a scale of hundreds only in the free 
shear layers, and clearly originating near the walls. In contrast, the solution with excessive ambient v has a high 
plateau away from the airfoil. The S-A destruction term activates in response to the high plateau level and strongly 
reduces v/v, from 5000 gradually towards 0 near the wall. It is still excessive, of course, and the eddy-viscosity fields 
are grossly different, particularly in the free shear layers. Here we discuss the reasons why the excessively high v 
degrades the solution. The approximate destruction equation is: u • V(l/F) = 2c w \/d 2 , where d is the distance to the 
wall (because in irrotational regions r 1, so that /,,, = c w 3 = 2). For a streamline directly approaching the airfoil 
without deceleration, its solution is 


1 1 ^ ^Cuji 

v vfs U d 


( 10 ) 


which limits v from above with Ud/2c w \ (recall that c w 1 = 3.2). This explains why the plateau at v = 5.5 x 10 ~ a UL 
starts dropping near d = 0.05, ahead of the slat, in Figure 4b. The drop is then faster than in Equation 10, because 
the local velocity u is rapidly dropping below U. For streamlines that flow past the airfoil, the velocity is closer to 
U, but the destruction effect is cumulative, and a kind of super-boundary layer forms. For a gap flow that would be 
entered by fluid without deceleration, the rough upper bound U d)2c w \ gives a gap Reynolds number no lower than 
/ lc ,„ 1 , or 12.8; typically it is well above 12.8 thanks to destruction accumulated along the airfoil. In figure 4b, the eddy 
viscosity in the flap gap does not exceed v/v = 750; compare with vfs/ v — 5000 and U d/2c w \/v k, 7000. The gap 
Reynolds number is about 180. With the normal inflow value in figure 4a, the gap Reynolds number is just 375 if based 
on the peak eddy viscosity, which propagates from the recirculation region in the flap cove. These two gap values for 
the different ambient levels are not hugely different, showing how the S-A model corrected the excessive inflow value 
to a good extent, for the gap. It also did rather well for the stagnation point. These are favorable consequences of a 
term which was calibrated for destruction in the logarithmic layer. However, the destruction term offers the free shear 
layers little protection: the peak eddy viscosity in the slat wake is artificially raised from about 2000 to 3000 when 
freestream vfs/ v goes from 5 to 5000. 

Figure 5 illustrates the end effect of the variations on the lift coefficient. Again, v/v needs to reach 500, or about 
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Figure 3. Multi-element airfoil flow fields using S-A, M=0.2, Re c = 9 million, showing (a) Mach contours near the flap gap for v /v = 5, 
(b) same for u j y = 50000, and (c) profiles of Mach number along x/c = 0.88, for u/v = 5 to 5000 in factors of 10. 



Figure 4. Multi-element eddy-viscosity fields ut/v using S-A, M=0.2, Re, = 9 million, showing (a) inflow vpsl v = 5; (b) same for 
vfs! v = 5000. Courtesy S. Allmaras. 
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Re c / 2 x 10 4 , before the effect on lift is appreciable (3% with v/v = 500, and 12% with v/v = 5000). The estimate 
given earlier was about Re c / 10 4 , therefore quite close. 



Figure 5. Lift coefficient of McDonnell-Douglas three-element airfoil with a range of values for ambient eddy viscosity. Chord Reynolds 
number 9 x 10 6 . Courtesy S. Allmaras. 

Turning our attention to the leading edge, we find that the values of 5 /L are quite Reynolds-number dependent 
since the laminar thickness is of order \J vr/U with r the leading-edge radius, and the turbulent thickness is often 
not much larger. However, the constraint on e or lo will be more severe if based on the thicker regions, where 5 is 
a few percent of chord. Free-shear regions such as the slat wake have similar if not larger thicknesses, and smaller 
velocity differences, which weakens their resistance. The merger of the slat wake with the main-element wake has 
much influence on the flow at some angles of attack, and therefore this wake cannot be “sacrificed.” This is not to 
imply that any of the present models approach perfection in predicting this merger; only that the solutions should be 
consistent between different codes and users. 

C. Recommended Values 

The considerations above lead to going-in upper bounds for the ambient primary-variable values k ss 1 x 10~ 4 U 2 , 
e ss 4.5 x 10 _5 C/ 3 /c, lo « 5 U/c. Here, c = L is the airfoil chord. Use of these upper bounds implies v t ss 2x 10 ~ 5 Uc, 
which is small enough for flap gaps and well-developed boundary layers, but possibly too large for the leading-edge 
boundary layers. This may not strongly influence the rest of the flow, because the acceleration weakens the memory of 
the boundary layer. Still, some applications such as laminar-flow control and icing are dependent on high accuracy in 
the leading-edge region. This suggests lowering the allowable k, to say, k ~ lx 10 ~ 6 [/ 2 , while preserving lo « 517 /c. 
This gives an upper bound on e of sa 4.5 x 10 ‘U 3 /c, and u t ss 2 x 10 7 Uc. 

Note that this ambient value of lo (or e/C^k) is unreachable if free decay takes place. An accurate solution for free 
decay sets the value of k/e in the neighborhood of 100c/17 for a 50c domain, or lo « 0.1 U/c. Such a low value is not 
required by comparison with the boundary layer (fig. 2), and is undesirable for the following reason: the SST model 
can encounter a particular problem when the ambient value of v t /v is large and lo is set too small. The SST limiter 
chooses the maximum of ciiu> and f 1F 2 (1> the vorticity). Outside the boundary layer, Cl — > 0 and F 2 — » 0, and it is the 
intention that the model choose a\LO in this region. However, if lo is too small and residual vorticity is present, the two 
terms compete and the result can be an irregular eddy-viscosity distribution, as shown in Fig. 6. The main problem 
with this erratic behavior is that convergence can stall. Some implementations of the SST model use the strain rate 
instead of 0, which is neutral in thin shear flows, but this is not the reference formula. In our tests with the NACA 
0012 airfoil, values of v t /v > \ and loc/U < 3 or so caused such problems. Imposing floor values are quite helpful 
in dealing with this problem. Figure 3 suggests that values of loc/U up to 100 would be more than safe, so that the 
useable range is not too narrow. 

To summarize, the following ambient freestream turbulence values (in the region near the body) are recommended 
for general use for most typical external aerodynamic applications with two-equation turbulence models: k/U 2 = 
1 x 10^ 6 (or Tu near « 0.1%), ec/U 3 = 4.5 x 10~ ' , loc/U = 5. These values yield an ambient eddy viscosity ratio 
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Figure 6. Example of poor solution behavior of SST model outside the boundary layer edge, due to too-low ambient uiL/U = 0.2, 
vt/v = 12.5, NACA 0012 airfoil, a = 0°, Re c = 10 7 . 


of v t /v ss 2 x 10“ ' x Re c . Thus, for example, for Re c = 10 million, v t /v « 2; and for Re c = 1 million, v t /v ss 0.2. 
Again, these ambient levels are not achievable if the turbulence variables are allowed to decay naturally according to 
the turbulence equations. One easy way to achieve these ambient levels in practice is to impose them as both inflow 
and floor values. 

For S-A, the values v/v = 3 and 5 are well-placed for fully-turbulent behavior, with 3 somewhat preferable at 
the lowest Reynolds numbers, since it gives v t /v ss 0.2, well below 1. Higher values may slightly help convergence 
at high Reynolds numbers by smoothing out the edges of turbulent regions, but not enough to motivate giving up the 
simple recommendation. 

Figure 7 explores the S-A flow field for a multi-element configuration using vps/ v = 5, with emphasis on the 
stagnation regions. An exponential scale is used in this figure for eddy viscosity. It is quite smooth on this grid, 
although with a little noise in the 5-to-10 range. Although not shown, a computation on a coarser grid (about 1/4 
the number of points) was not smooth for the low values of v up to about 50, and it also had significant negative 
excursions of v in the flap gap, traced to irregularities in the adaptive unstructured grid. The fine grid has negative 
values only very far downstream (and in this code, when v < 0, vt is set to 0). Multi-block structured grids generate 
their own disruptions, but they appear to be milder, and confined to block comers and small regions with highly 
stretched and distorted cells. There is no evidence here that the imperfect solution is harming the “physics” of the 
model, but a blanket statement can certainly not be made. For instance, contact between unphysical “hot spots” and 
the flap upper-surface boundary layer could precipitate transition in a grid-dependent manner. 

Also striking in this computation is the fact that the model does not depart from laminar for quite some distance in 
some of the boundary layers. For instance, the lower stagnation point on the slat is near x = —0.035, and appreciable 
eddy viscosity does not grow until the adverse pressure gradient at x = —0.085. This supports an earlier comment 
about the boundary-layer thicknesses in laminar and turbulent flow in stagnation regions being close (even at a high 
chord Reynolds number). The boundary layer is vigorously thinned out by the acceleration, yet the model remains 
poised to grow turbulent values, because the ambient value was chosen for that. The flap lower-surface boundary layer, 
similarly, is free of significant turbulence within the frame shown; it grows very gradually, to reach 27 at the trailing 
edge. The Turbulence Index i t (ref. 4) is about 0.5 at the stagnation points, but rises to 1 rapidly. 

Figure 8 is a similar figure for the SST model; the left two frames can be directly compared with those in Figure 
7. The three turbulence variables are shown in exponential scales again, with the darkest blue for the ambient values. 
The recommended ambient values (also used as floor values) appear very effective, in that they are smaller than the 
boundary-layer values almost everywhere along the wall, yet transition is occuring reliably. Note mild differences 
between the S-A and SST models in the slat cove: the peak eddy-viscosity values are 1650 and 1800, respectively, and 
the SST recirculation bubble is accordingly wider. Conversely, the slat-upper-surface boundary layer near x = —0.025, 
S-A peaks at 314 and SST at 277 , and the S-A boundary layer is marginally thicker. Similar runs with the chord 
Reynolds number ranging from 10 6 to 10 s were all successful. 
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Figure 7. Detail of S-A solution, McDonnell-Douglas three-element airfoil. Upper frames, eddy viscosity; lower frames, velocity magnitude. 
Chord Reynolds number 9 x 10®, a = 24°, upg/v = 5. Unstructured grid, courtesy S. Allmaras. 
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Figure 8. Detail of SST solution, McDonnell-Douglas three-element airfoil. Upper left, eddy viscosity vt/v\ upper right, kps/U 2 ; lower 
left, velocity magnitude; lower right, ujc/U. Re c = 9 x 10 6 , a = 19°, kps/U 2 = 1 x 10 -6 , ujpsc/U = 5, with floor values. 
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Finally, Fig. 9 demonstrates how use of the recommended ambient values in the SST model can reduce the grid- 
dependence of apparent transition location for an airfoil computation. In Rumsey, 1 the grid dependence was first 
demonstrated for typically-used freestream turbulence levels, but a satisfactory explanation and cure were not found. 
Fig. 9(a) shows the sort of grid-dependent transition seen. Although relatively small, this grid-dependence could 
conceivably impact conclusions drawn from grid sensitivity studies, for example. In Fig. 9(b) we find that by applying 
the currently recommended freestream levels, and limiting via floor values (at the same levels), the grid-dependence 
in transition location is largely eliminated. 
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Figure 9. Skin friction coefficient near leading edge of NACA 0012 airfoil at Re c = 10 7 , with SST model; (a) grid-dependent solution from 
Rumsey 1 using typical freestream turbulence levels and no floor values; (b) solution using recommended freestream turbulence levels and 
floor values. 


IV. Summary 

Inflow conditions for CFD with simple turbulence models in large domains were studied from theoretical and 
pragmatic points of view, with the following outcomes: (1) simplistic and erroneous conceptions in the community 
were discussed; (2) rates of free decay for widely-used two-equation model levels were shown to be under-resolved 
on typical grids; (3) excessively high values of eddy viscosity were demonstrated to contaminate the flow field in 
non-turbulent regions; (4) the use of floor values, which greatly facilitate the control of the ambient levels in two- 
equation models, was advocated; (5) rational lower and upper bounds for these ambient values were provided; and 
(6) specific ambient turbulence values were recommended as safe starting points, located within rather wide intervals 
of low sensitivity. The practice of floor values, which conflicts with the two fundamental turbulence equations, was 
supported by theoretical considerations on the relevant range of turbulent motions. Thanks to having a single equation, 
the Spalart-Allmaras model readily provides an equivalent level of control, plus the ability to steer the boundary layers 
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to be laminar or turbulent, without any special measures. If accepted, the present recommendations will spread in CFD 
codes, lowering the user’s burden and re-work, improving convergence and consistency between different users, and 
improving the accuracy of CFD solutions by removing unseen contaminations and helping users employ the models 
precisely as their authors intended. 
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